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Abstract: Adequate modeling of light propagation in the head is important 
to predict the sensitivity of MRS signal and the spatial sensitivity profile of 
source-detector pairs. The 3D realistic head models of which the geometry 
is based upon the anatomical images acquired by magnetic resonance 
imaging and x-ray computed tomography are constructed to investigate the 
influence of the frontal sinus on the MRS signal and spatial sensitivity. 
Light propagation in the head is strongly affected by the presence of the 
frontal sinus. The light tends to propagate around the frontal sinus. The 
influence of the frontal sinus on the sensitivity of the MRS signal to the 
brain activation is not consistent and depends on the depth of the frontal 
sinus, the optical properties of the superficial tissues and the relative 
position between the source-detector pair and the frontal sinus. The frontal 
sinus located in the shallow region of the skull tends to reduce the 
sensitivity of the MRS signal while the deep frontal sinus can increase the 
sensitivity of the MRS signal. 
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1. Introduction 

Near-infrared spectroscopy (NIRS) and optical imaging have been used to non-invasively 
measure cerebral oxygenation and hemodynamic changes associated with brain activation 
[1,2]. They have been widely applied to a variety of brain function studies [3-7] and clinical 
measurements [8-13]. In NIRS measurements, the source and detector probes are attached to 
the scalp to detect the change in optical density caused by the hemodynamic change in the 
brain. Consequently, the light needs to pass through the superficial tissues such as, scalp, skull 
both before and after passing through the brain tissue. It is important to obtain the sensitivity 
of the NIRS signal to the absorption change in the brain and the volume of tissue sampled 
with a particular source-detector pair, however, light is strongly scattered by the tissues, thus 
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the path of light propagation in the head is ambiguous. Since light intensity in the head cannot 
be non-invasively measured by experiments, light propagation in the head has been predicted 
by simulations. 

Light propagation in the head is affected by the heterogeneity in optical properties of the 
tissues. Especially, cerebrospinal fluid (CSF) that has both low scattering and absorption 
coefficients has a strong effect on light propagation in the brain [14-16]. The influence of a 
CSF layer cannot be ignored in any modeling of light propagation in the head. The boundaries 
between various tissues in the head have complex geometry. The brain surface is folded with 
sulci filled with the CSF. The geometry of the tissue boundaries is also important for adequate 
modeling of light propagation in the head. Recently, the realistic head models, the structures 
of which are based upon magnetic resonance (MR) head images has been used to predict light 
propagation in the head for NIRS and optical imaging [17-20]. These head models basically 
consisted of five types of tissues that imitate the scalp, skull, CSF, gray matter and white 
matter. A more sophisticated model, for instance, which includes the extra-cerebral 
vasculature, has been constructed to investigate the influence of the large pial vessels, such as 
sagittal sinus, on the diffuse optical imaging [21]. 

The frontal sinus (FS) in the skull is an optically void region. The analysis of light 
propagation in a simplified head model with a void region mimicking the FS indicated that 
light propagation in the brain is affected by the FS [22]. However, the FS tends to be ignored 
in the realistic head models based on the MR images because there is very little contrast 
between the FS and skull regions in the MR image to obtain precise geometry of the FS. 

In this study, the 3D realistic head model including a FS is constructed in order to 
investigate the influence of a FS on the NIRS measurements. The FS can be clearly observed 
in images acquired by x-ray computed tomography (CT). The anatomical geometry of the 
head model is based upon MR and x-ray CT images. Three additional models including the 
artificially modified FS are also constructed and light propagation in these head models is 
compared to discuss the influence of the FS on the NIRS signal and the spatial sensitivity 
profile. 

2. Method of simulation 

2.1. Construction of the realistic head models 

The geometry of the realistic head models is based on the anatomical images obtained by MR 
imaging and x-ray CT. The original 3D MR images contain 166 planes with 1.0 mm intervals 
and each plane consists of 256 x 256 pixels with square sides of 1.0 mm. The original 3D CT 
images contain 381 planes with 0.5 mm intervals and each plane consists of 512 x 512 pixels 
with square sides of 0.41 mm. Figures 1(a) and 1(b) show the Tl -weighted MR image and CT 
image in the same plane of a slice, respectively. In the Tl -weighted MR images, the contrast 
among the grey matter, white matter and CSF is sufficient enough to identify these 
boundaries. It is difficult to find the boundaries between the skull and the adjacent regions 
because bone and water (CSF) produce a low signal while the bone marrow and fat in the 
scalp produce a high signal. Since the air produces a low signal in the Tl -weighted MR 
images, the region of the FS cannot be recognized in the images. In the CT images, the 
boundaries between the skull and adjacent regions are easily distinguishable and the FS is 
clearly observed in the frontal bone. 

The head models consist of 6 regions, which include the air, scalp, skull, CSF 
(subarachnoid space), gray matter and white matter. Since the spatial resolution of the MR 
and CT images are different, the CT images were realigned by SPM8 software [23] to 
coregistrate images of the MR and CT. The MR and CT images were binarized by the 
conventional region growing algorithm [24] for the 3D grayscale images. The outer surface of 
the whole head, gray matter and white matter were identified from the MR images while the 
skull and air in the FS were identified from the CT images to determine the boundaries of 
each type of tissues. An axial slice of the head model segmented from the MR and CT images 
and the 3D segmented geometry of the head model are shown in Figs. 2(a) and 2(b), 
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(a2) 




(bl) 



Fig. 1. (al) An axial slice of MR image. (a2) A sagittal slice of MR image, (bl) An axial slice 
of CT image. (b2) A sagittal slice of CT image. The MR and CT images are in the same plane. 
White lines in (al) and (a2) indicate the plane of the sagittal and axial slices, respectively. 




Fig. 2. (a) An axial slice of the segmented image of the realistic head model, (b) 3D segmented 
geometry of the realistic head model. The outer surfaces of the whole head, FS and brain are 
only represented. 



respectively. The outer surfaces of the whole head, FS and brain are represented in the 3D 
segmented geometry. The head models consist of 256 x 256 x 106 cubic voxels with side 1 
mm. Each voxel is specified with its optical properties according to the results of the tissue 
segmentation. The optical properties of each type of tissues in the head models at the 830 nm 
wavelength were chosen from the reported data [16,25-27] (parameter set (I)). Actually, a 
wide range of different values in the optical properties has been reported for the same types of 
tissues. Light propagation in the head models with two additional parameter sets was 
predicted in order to investigate the dependence of the influence of the FS on the optical 
properties of the head. Extreme values were chosen based on the literature [21,28,29] as 
parameter set (II). The scattering coefficient of all the types of tissues in parameter set (II) is 
lower than those in parameter set (I). In parameter set (III), the optical properties of the scalp 
and skull are the same as those in parameter set (II) while the optical properties of the other 
parts of the head are the same as those in parameter set (I). The optical properties of each type 
of tissues in the head models are listed in Table 1 . 
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Table 1. Optical properties of each type of tissues in the head model 





Parameter 
set 


Scalp 


Skull 


FS 


CSF 


Gray 
matter 


White 
matter 


Transport scattering 


I 


1.84 


1.48 


0.01 


0.24 


2.12 


8.82 


coefficient (mm -1 ) 


II 


0.66 


0.86 


0.01 


0.01 


0.7 


0.7 


III 


0.66 


0.86 


0.01 


0.24 


2.12 


8.82 


Absorption coefficient 


I 


0.025 


0.025 


0.001 


0.004 


0.039 


0.017 


(mm" 1 ) 


II 


0.019 


0.014 


0.001 


0.003 


0.020 


0.080 


III 


0.019 


0.014 


0.001 


0.004 


0.039 


0.017 



2.2. Modification of the FS 



Three additional models were constructed to investigate the influence of the FS on the 
sensitivity of the MRS signal. In the additional models, the geometry of the FS was only 
modified from that in the realistic head model shown in Fig. 2 (the model with the whole FS). 
Figure 3(a) shows an axial slice of the frontal area of the model with the whole FS. In the 
model without the FS, the FS is replaced by the skull as shown in Fig. 3(b). In the model with 
the shallow FS, half of the volume in the deeper region of the FS was replaced by the skull 
and that in the shallow region was replaced by the skull in the model with the deep FS as 
shown in Figs. 3(c) and 3(d), respectively. 

(a) 
(c) 

Fig. 3. An axial slice of the head models for evaluation of the influence of the FS on light 
propagation in the head, (a) The model with the whole FS, (b) the model without the FS, (c) the 
model with the shallow FS, (d) the model with the deep FS. 

2.3. Analysis of light propagation in the head models 

Light propagation in the head models was calculated by use of a Monte Carlo method [30- 
32]. Refractive index mismatch was considered only at the interface between the scalp and air. 
The refractive indices of air and tissue were 1.0 and 1.4, respectively. A fraction of the photon 
weight absorbed at each interaction site was accumulated in each voxel in time-gates of width 
50 ps to calculate the time-resolved photon fluence. When a photon packet reached the 
detector position, its ultimate weight was accumulated to calculate the detected intensity. The 
path length of the detected photon packet in each voxel and each type of tissues weighted by 
the ultimate weight was accumulated to calculate the spatial sensitivity profile and partial 
optical path length in each type of tissues. The partial optical path length can be used as an 
index of the sensitivity of the MRS signal to global absorption change in a particular type of 
tissues while the spatial sensitivity profile indicates the volume of tissue sampled by the 
source-detector pair of the NIRS instrument. The positions of the sources and detector are 
shown in Fig. 2(b). The distance between each source-detector pair is 30 mm. 
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3. Results 



3.1. Influence of the FS on light propagation in the head models 

The time-resolved photon fluence maps in the head models with the whole FS and without the 
FS shown in Figs. 3(a) and 3(b), respectively, were predicted in order to evaluate the 
influence of the FS on light propagation in the head. A pulse of light was irradiated at source 
position A in Fig. 2(b) and the photon fluence maps during every 50 ps were recorded up to 
2.5 ns. A frame from a movie (Media 1) of the time resolved photon fluence map in the head 
models with parameter set (I) is shown in Fig. 4. When the photon fluence reaches the FS at 
about 250 ps, the photon fluence around the FS rapidly increases. This causes the distortion of 
the photon fluence map and the distribution of the photon fluence is broadened in the direction 
of the FS. The penetration of the photon fluence tends to be confined to the gray matter and 
the photon fluence laterally spreads in the gray matter in both models. This phenomenon is 
caused by the low- scattering and low-absorption characteristics of the CSF. 











400 ps 















Fig. 4. Frame from a movie (Media 1) of light propagation in the head model with the whole 
FS and without the FS. The red region indicates the FS. 

3.2. Influence of the FS on the NIRS signal and spatial sensitivity profile 

Light propagation in the four head models with and without the FS was predicted in order to 
evaluate the influence of the FS on the NIRS signal. The source at position A or B and 
detector were placed on the forehead of the model as shown in Fig. 2(b). In the case of source 
A, the FS existed inside the distribution of the light path connecting the source and detector, 
while the FS was almost outside of the light path in the case of source B. 

The spatial sensitivity profiles for the source at position A and detector in the four head 
models with parameter set (I) are shown in Fig. 5. The spatial sensitivity profile is distorted by 
the FS and the distortion depends on the depth of the FS. In the models with the FS, the spatial 
sensitivity in the scalp around the midpoint between the source and detector is significantly 
reduced and the sensitivity around the FS increased. The depth profile of the spatial sensitivity 
along the line in Fig. 5 is shown in Fig. 6. The whole FS and shallow FS reduce the 
penetration depth of the spatial sensitivity profile. The deep FS increases the penetration depth 
of the spatial sensitivity profile, however, the penetration depth in the gray matter is almost 
the same as that in the model without the FS. The spatial sensitivity profiles on the brain 



#170737 - $15.00 USD Received 18 Jun 2012; revised 29 Jul 2012; accepted 29 Jul 2012; published 14 Aug 2012 
(C) 2012 OSA 1 September 2012 / Vol. 3, No. 9 / BIOMEDICAL OPTICS EXPRESS 2126 





Fig. 5. The influence of the FS on the spatial sensitivity profile for the source at position A and 
detector, (a) The model with the whole FS, (b) the model without the FS, (c) the model with the 
shallow FS, (d) the model with the deep FS. The line in the figures represents the position of 
the depth profile shown in Fig. 6. 



M 

o3 



^ 0.08 - 



•B 0.06 - 



0.04 



•2 0.02 



CSF ; gray matter 




5 10 15 20 

depth from the surface of the head (mm) 

Fig. 6. The depth profile of the spatial sensitivity profile in the four head models. The position 
of the profile is shown in Fig. 5. 



surface in the four models with parameter set (I) are shown in Fig. 7. The dashed line 
indicates the position of the FS. A hypothetical position of the FS is also shown in the model 
without the FS for reference. The spatial sensitivity profile on the brain in the model without 
the FS shown in Fig. 7(b) is almost symmetrically distributed with respect to the line 
connecting the source and detector while that in the model with the whole FS shown in Fig. 
7(a) is slightly shifted toward the lower side and the peak of the sensitivity is reduced. The 
spatial sensitivity is narrowed in the model with the shallow FS while the spatial sensitivity is 
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(c) (d) 



Fig. 7. The influence of the frontal sinus on the spatial sensitivity profile on the brain surface, 
(a) The model with the whole FS, (b) the model without the FS, (c) the model with the shallow 
FS, (d) the model with the deep FS. Dashed line indicates the position of the FS. 



Table 2. Intensity and path length of detected light for the head models with different 

types of FS 



Type of 


Detected 


Total path 






Partial optical path length (mm) 




frontal sinus 


intensity 


length (mm) 


Scalp 


Skull 


FS CSF Gray matter 


White matter 


(a) Whole FS 


3.4 x l(T 8 


166 


88.4 


67.7 


37.6 7.81 2.06 


0.253 


(b) Without FS 


7.5 x 10" 9 


230 


115 


101 


0 11.3 2.71 


0.397 


(c) Shallow FS 


2.7 x 1(T 8 


179 


91.6 


80.7 


29.0 5.23 1.35 


0.163 


(d) Deep FS 


1.3 x l(T 8 


210 


100 


94.7 


16.6 12.4 3.01 


0.404 



broadened toward the lower side in the model with the deep FS. These results indicate that the 
increase in photon fluence around the FS affects the spatial sensitivity profile of the NIRS 
measurements. 

The tendency of the influence of the FS can be quantitatively evaluated by the intensity 
and optical path length of the detected light shown in Table 2. The partial optical path length 
in the scalp is reduced while that in the skull including the FS is increased in all the models 
with the FS. This indicates that the light which propagates through the FS tends to be detected 
by the detector. Since light propagates in straight line without absorption in the FS, the 
detected intensity is increased and the total optical path length is decreased by the presence of 
the FS. It should be noted that the partial optical path length in the FS is unexpectedly long 
although no scattering can occur in the FS. This indicates that many photons scattered in the 
skull passes through the FS without absorption. These photons increase the photon fluence 
around the FS. The influence of the FS on the partial optical path length in the gray matter 
depends on the depth of the FS. The deep FS increases the partial optical path length in the 
gray matter while the shallow FS decreases the partial optical path length in the gray matter. 

Figure 8 shows the ratio of the partial optical path length in the gray matter for the head 
models with the FS to that for the head model without the FS. The results for the models with 
the three different parameter sets are compared in order to investigate the dependence of the 
influence of the FS on the optical properties of the head. In the head models with all the 
parameter sets, the deep FS increases and the shallow FS decreases the partial optical path 
length in the gray matter. However, in the case of the whole FS, the partial optical path length 
in the gray matter for the models with parameter sets (II) and (III) are increased by the 
influence of the FS while that in the model with parameter set (I) is decreased. The tendency 
of the influence of the FS in the head model with parameter set (III) is almost the same as that 
in the head model with parameter set (II). This result indicates that the influence of the FS on 
the NIRS signal depends on not only the depth of the FS but also the optical properties of the 
superficial tissues such as the scalp and skull. 
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Fig. 8. The influence of the optical properties of the head model on the partial optical path 
length in the gray matter. The partial optical path length in the gray matter for the head models 
with the FS is normalized by that for the head model without the FS. Each parameter set is 
listed in Table 1 . 

Figures 9(a) and 9(b) show the spatial sensitivity profiles for the source at position B and 
detector in the models (parameter set (I)) with the FS shown in Fig. 3(a) and without the FS 
shown in Fig. 3(b), respectively. The spatial sensitivity profile in the skull is significantly 
distorted toward the FS as shown in Fig. 9(a). Consequently, the spatial sensitivity profile on 
the brain surface is broadened toward the direction of the FS. The influence of the FS on the 
sensitivity of the MRS signal depends on the relative position of the source-detector pair and 
the FS. 





Fig. 9. The influence of the FS on the spatial sensitivity profile for the source at position B and 
detector, (a) The model with the whole FS. The FS is almost outside of the distribution of the 
light path connecting the source and detector, (b) The model without the FS. 



4. Discussion 



Light propagation in tissue is known to be strongly affected by the presence of voids and low- 
scattering regions. The light tends to propagate around the void region, and consequently, the 
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spatial sensitivity around the void region is increased. In the brain activation measurements by 
MRS, the low-scattering CSF affects light propagation in the head and increases the partial 
optical path length, i.e., the sensitivity of the MRS signal, of the gray matter neighboring the 
CSF [14—16]. The FS also increases the spatial sensitivity around the FS, and consequently, 
decreases that apart from the FS. Since the brain is not next to the FS, the influence of the FS 
on the sensitivity of the MRS signal is not consistent and is more complicated than that of the 
CSF. The influence of the FS on the sensitivity of the MRS signal depends on the depth of the 
FS and the optical properties of the scalp and skull. The shallow FS decreases the sensitivity 
of the MRS signal while the deep FS increases the sensitivity. The depth at which this 
inversion based on the influence of the FS occurred is changed by the optical properties of the 
scalp and skull. The optical properties of the tissues depend on the wavelength and the wide 
range of different values for the same type of tissues is obtained from the literature. The 
influence of the FS also depends on the relative position of the source-detector pair and FS. 
When the FS exists almost outside of the distribution of the light path connecting the source 
and the detector, the spatial sensitivity profile tends to be broadened. Furthermore, the 
realistic head models in this study do not consider the presence of the extra-cerebral 
vasculature, such as sagittal sinus. The influence of the sagittal sinus on the MRS signal has 
been reported [21] and the sagittal sinus is located behind the middle of the FS. Careful 
consideration is needed to estimate the influence of the FS on the sensitivity of the NIRS 
signal. 

The FS tends to locate below the eyebrows, and in this case, the influence of the FS on the 
NIRS measurements of the brain activation in the frontal area may be minimized by the 
optimal probe arrangement. However, the depth and size of the FS is different between 
individuals. It is necessary to analyze light propagation in the 3D head model, the anatomical 
structure of which is based on the structural image of the individual if you need to obtain 
details of the influence of the FS on the MRS measurements. 

5. Conclusions 

The 3D realistic head models of which the structure is based on the MR and x-ray CT images 
are constructed to investigate the influence of a FS on the brain function measurements by 
near-infrared spectroscopy and imaging. Light propagation in the head model with the FS of 
which the geometry was based on the CT images was compared to that in the head models 
with the modified FS to investigate the influence of the FS on the MRS signal. Light 
propagation in the head is strongly affected by the presence of the FS. The distribution of the 
photon fluence is broadened in the direction of the FS, consequently, the spatial sensitivity 
around the FS tends to be increased. Since the influence of the FS on light propagation in the 
brain depends on the depth of the FS and the optical properties of the scalp and skull, the 
influence of the FS on the sensitivity of the MRS signal to the brain activation is not 
consistent across the source-detector pairs. The shallow FS tends to reduce the sensitivity of 
the MRS signal, however, the deep FS can increase the sensitivity of the MRS signal. 
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